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Abstract. - The Loschmidt echo and the purity are two quantities that can provide invaluable information 
about the evolution of a quantum system. While the Loschmidt echo characterizes instability and sensitivity 
to perturbations, purity measures the loss of coherence produced by an environment coupled to the system. 
For classically chaotic systems both quantities display a number of ~ supposedly universal ~ regimes that 
can lead on to think of them as equivalent quantities. We study the decay of the Loschmidt echo and the 
purity for systems with finite dimensional Hilbert space and present numerical evidence of some fundamental 
differences between them. 



Some of the latest breakthroughs in theoretical and exper- 
imental quantum physics have permitted among other things 
to explore and manipulate new states of matter - like Bose- 
Einstein condensates - and also manipulate small numbers of 
atoms or ions making it possible to test some of the assertions 
of relatively new areas of research like quantum information. 
Two of the main problems that affect the achievement of such 
advances are uncontrolled coupling to an environment and ir- 
reversibility caused by sensitivity to small perturbations of the 
quantum evolution. 

The presence of a coupling to an external bath introduces 
decoherence (T\. By definition decoherence washes out inter- 
ference terms due to quantum superposition. One way of char- 
acterizing the decrease of the interference terms caused by de- 
coherence is by measuring the purity of the system as a function 
of time. For classically chaotic systems it was conjectured ^ 
and numerically shown 123 that for a certain range of val- 
ues, the exponential decay of purity is independent of the cou- 
pling strength and is characterized by the Lyapunov exponent 
of the classical counterpart. Complementarity, to characterize 
irreversibility and instability arising from the chaotic nature of 
systems the Loschmidt echo (LE) - also known as fidelity - has 
been used ||5H9l. The idea is to study the overlap as a function 
of time of two states evolving with slightly different evolution 
operators characterized by some perturbation parameter E. To 
avoid the singularity of a particular initial condition the LE is 
computed by averaging over many initial states. The averaging 



can be treated in analogy to the effects of decoherence and it 
was claimed flOl - or expected ifTTl - that at least for classi- 
cally chaotic systems, since they exhibit the same decay rates, 
they provide essentially the same information. Both the purity 
and the LE are of capital importance in experiments because 
they are readily measurable. In quantum information fidelity 
helps determine the accuracy of quantum gates and state prepa- 
ration. Purity on the other hand characterizes how much the 
system is coupled to an environment or how much two parties 
of a system are entangled with each other. 

In the present contribution the aforementioned quantities 
are explored for systems with finite dimensional Hilbert 
space, and which have a classically chaotic counterpart. We 
present numerical evidence that, contrary to previous results 14] 
ilOl, significant differences can exist between the -supposedly 
universal- behaviors of the LE and the purity. We will focus on 
the asymptotic regime as a function of perturbation [and deco- 
herence] strength, both in the small perturbation regime as well 
as for larger perturbations. Differences in the short time regime 
were already studied elsewhere [see e.g. 1 12]]. While for small 
perturbations the LE shows the expected quadratic regime, we 
show that the decay rate of the purity, as a function of the cou- 
pling strength, depends -strongly- on the type of environment 
affecting the system. Furthermore, in the strong perturbation 
regime, the LE presents an oscillatory behavior that can mask 
the Lyapunov decay [ 13 - 16 1. A measurement of the fidelity de- 
cay in these regimes can thus give results that are far from the 
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expected. Besides, the Lyapunov decay for the purity is only 
observed for special types of environment. 

The systems we consider are quantum maps on the torus. 
Quantum maps provide a useful tool to understand universal 
properties of quantum chaotic systems - e.g universal spectral 
statistics. In addition, there exist efficient quantum algorithms 
for some quantum maps that can be implemented with a small 
number of qubits ITtHTsI . This makes them ideal testbeds 
for current quantum computers. The quantized torus has as- 
sociated an dimensional Hilbert space with Planck constant 
H = \/2nN, and the position basis {^i}^^ and momentum ba- 
sis {pi}^^ ' are related by the discrete Fourier transform (DFT). 
We consider quantum maps U whose classical counterpart are 
chaotic. For simplicity, we use maps whose evolution operator 
for one iteration can be written as follows 



J2nNT{p)-a7zNV{q) 



(1) 



These kind of maps can be efficiently computed using the DFT 
(through the fast Fourier transform) and some have been imple- 
mented experimentally - e.g [19J- and have efficient quantum 
algorithms - e.g. II17II18I . The corresponding classical map is 
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In particular for numerical calculations we use the cat map per- 
turbed with a non-linear shear 



q' 



p + aq + 27zk(cos[2nq] — cos[4;r^]) 
q + bp' 



(modi), (3) 



with a,b integers. This map is chaotic with largest Lyapunov 
exponent A w ln((2 + ab + ^Jab{A + ab))/2)/2, for ;t < 1. 
The perturbation destroys the symmetries related to the arith- 
metic nature of cat maps - which are responsible for non- 
generic spectral statistics. 

Let us first consider the study of the LE. For pure states the 
LE is defined as 



M{t)^\{w\ullui\w)?. 
We define the perturbation strength as 
!. = !<! -k. 



(4) 



(5) 



In what follows t is an integer that represents the number of 
times that the map U is applied. Eq. (|4]l is dubbed "echo" be- 
cause it measures the overlap between a state evolved forwards 
up to time t with U and then backwards with the slightly per- 
turbed operator t/j;. It can also be seen as a measure of the 
separation of two, initially identical states, evolved forwards 
with two slightly different evolution operators. If the classical 
dynamics is chaotic, there are three well identified regimes for 
the LE as a function of time: parabolic or Gaussian for very 
short times; exponential for intermediate times followed by a 
saturation depending on the effective Hilbert space size. Here 
we focus on the decay rate Fle as a function of E for the ex- 
ponential decay regime. The decay rate can be extracted from 



the smooth curves [see Fig. [T] obtained after performing an 
average over 1024 uniformly distributed intial coherent states, 
chosen randomly. 

Fig- m [left] shows the decay of the LE as a function of 
discrete time t for various values of E. We see that after a 
few steps the decay is exponential. As expected for small 
perturbation strength (e.g. E/fi = 0.65884, L31768, L6471) 
the decay rate is much smaller than A, but contrary to pre- 
dictions l^/Tl as we increase E (e.g l./h = 2.63536,2.96478) 
decay rates can reach values much larger than X. The Lya- 
punov decay appears for greater values of the perturbation (e.g. 
E/Zi = 6.5884). This complex behavior of the decay rate of 
the LE as a function of the perturbation is shown in detail in 
Fig. [U [right] where we plot the decay rate Fle as a function of 
the rescaled strength of the perturbation E/Zi. We tested results 
for two different versions of the map of Eq. Q with differ- 
ent Lyapunov exponent [(0) a = b ~2, X ^ ln(3 + 2\/2); (•) 
a = b = A, X = ln(9 + 4-\/5)]. We can see that for small pertur- 
bation strength the behavior is, as expected, Fle °^T? - usually 
called Fermi golden rule (FGR) regime. For larger perturbation 
strengths, the decay rate is not as commonly predicted in the 
literature [see [8 , 9| and references therein] -with some excep- 
tions, e.g. llT3HT6l - perturbation independent behavior We find 
oscillations behavior near the value X. These oscillations can 
be understood through the local density of states (LDOS). For 
finite dimensional Hilbert space the LDOS grows quadratically 
with the perturbation up to a point where it starts to oscillate. If 
the mean value of the oscillatory part is comparable or smaller 
than the classical Lyapunov exponent, then the oscillatory be- 
havior is reflected in the echo. If, on the contrary, the Lyapunov 
exponent is much smaller than the mean value of the oscilla- 
tions of the LDOS, then no oscillations are appreciated in the 
LE 1 16|. The important thing to remark is that, after the FGR 
behavior, the decay of the LE is not perturbation independent. 
This can explain the difficulty to find the Lyapunov regime in 
echo experiments 1 15 ]. 

We now consider the evolution of our system in the pres- 
ence of an environment. We explore the behavior of the pu- 
rity for different types of environments. Interaction between 
system and environment produces global state which is non- 
separable, i.e. entangled. Once we trace out the environment 
degrees of freedom the reduced density matrix obtained evolves 
non-unitarily with a consequent loss of coherence. One way to 
measure the effect of the decoherence produced by the envi- 
ronment is through the purity [see e.g. | |20) ] as a function of 
time 



^(0=tr(p2), 



(6) 



were p, is the reduced density matrix of the system. The pu- 
rity is basis independent and measures the relative weight of 
the non-diagonal matrix elements. It can be used to measure 
how entangled are two systems coupled together. If P(f ) = 1, it 
means that the global system can be factorized into two sepa- 
rate systems and there is no entanglement. On the contrary for 
maximally entangled states the reduced density matrix has min- 
imum purity and the state is maximally mixed. In the case of 
an dimensional system P{t) = 1 /N for a maximally mixed 
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Fig. 1: Left panel: The LE as a function of time for the map of Eq. {3} with a = b = A, N = (corresponding to • symbols on 
the right panel), and averaged over 2^'^ initial states uniformly distributed. The diferent values of the reseated perturbtion are Z/fi = 
(□)0.65884, (■)!. 31768, (A) 1.6471, (0)2.63536, (♦)2.96478, (a)6.5884. The slope of the dashed line is the Lyapunov A = ln[9 + 4V^] 
2.88727. The horizontal line is at the saturation In(A'). Right panel: Decay rate Fle of the LE as a function of the reseated strength of the 
perturbation Z/fi. The map is the quantum version of the perturbed cat [Eq. ((S}] with (0) a = b = 2; (•) a = b = 4. Other parameters are: 
k = 0.0002, A' = 2^", and 1024 randomly chosen initial states. The lines are: (dashed) (L/h)^. The horizontal dotted lines correspond to the 
Lyapunov exponents of the corresponding maps: (below) A = ln[3 + 2\/2] ~ 1.76275; (above) ln[9 + 4\/5] ~ 2.88727. The inset shows the 
same in log-log scale where the quadratic small-Z regime is best appreciated. 



state. As a function of time, after an initial short transient, 
the purity decays exponentially. Like the LE for long times it 
saturates to a minimum value given by h. We focus on the ex- 
ponential decay and the dependence of the decay rate on the 
coupling parameter. 

Instead of studying the evolution of system plus environment 
and then tracing the environment out, we model directly the 
effect of the environment as a map of density matrices, or su- 
peroperator which, for Markovian environment and weak cou- 
pling, can be written in Kraus operator sum form 121]. The 
decoherence models we use can be expresed as a weighed sum 
of unitary operations. 



p' = D.(p) = 



E C£{q,p)T^pPT, 

p,q=0 



t 



(7) 



where T^p are the translation operators on the torus, /?) is a 
function of q and p and e characterizes the strength. The Kraus 
form implies complete positivity and the trace is preserved if 
\ — I. Furthermore, as T „ are unitary, the identity 



is preserved, i.e. the map is unital. Although position and 
momentum operators are not well defined in finite dimensional 
Hilbert space, translations can be defined as cyclic shifts [22,]. 
In Ref. Il23l it is shown that a variety of noise superoperators 
can be implemented in the form of Eq. The interpretation 
is simple: with probability C£{q,p) every possible translation 
in phase space is applied to p (incoherently). The decoherent 
effect of is evident: suppose we have a Schrodinger cat 
state that exhibits interference fringes in the Wigner function. 
Eq. ^ written for the Wigner function of p results 

W\Q,P)=Y,Csiq,p)WiQ-q,p^P). (8) 
9.P 



Then this incoherent sum of slightly displaced Wigner func- 
tions, washes out fast oscillating terms leaving only the classi- 
cal part. 

The complete map with decoherence takes place then in 
two steps, the unitary followed by the nonunitary part p' = 
D£{UpU'). This is an approximation that works exactly in 
some cases, e.g. a billiard that has elastic collisions on the walls 
and diffusion in the free evolution between collisions. 

To model diffusive decoherence we can define 



ceiq^p) = ^exp 



q 



2(i) 



(9) 



periodized to fit the torus boundary conditions. We will call 
this model Gaussian diffusion model (GDM). Eq. dHJ in the 
continuous limit is a convolution of the Wigner function with a 
kernel Cg{q,p). For the GDM this corresponds to the solution 
of the heat equation with diffusion constant given by {Ne/ln)^ 



In Fig. |2] [left] we show the behavior of the purity as a func- 
tion of time for map of Eq. (O with a = b ^2, N = 800 in the 
presence of GDM for different values of e. The exponential 
decay is clearly observed. Moreover, as £ increases (from right 
to left in the lines) the decay becomes independent of the value 
of e and is given by the classical value X. 

In continuos Hilbert space and in the presence of GDM type 
decoherence, the decay rate of the purity exhibits two different 
regimes as a function of the coupling parameter e. For small 
values of e it is equivalent to that of the LE as a function of E 
14, 10|, i.e. the decay rate depends quadratically for small £. 
Then, after a critical value it becomes independent of the envi- 
ronment and results F^ = A |i2jj4J. For large enough e values. 
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Fig. 2: Left panel: Purity as a function of time for for the map of Eq. l|3ll with a = = 2, = 2'" for various values of the decoherence parameter 
e. Increasing £ corresponds for lines from right to left (or bottom to top). The dashed stright line has a slope equal to A = ln[(3 + 2\/2)]. The 
dotted line shows the saturation at In(A'). Right panel: Decay rate (given by the slopes of the curves in the left panel) of the purity as 
a function decoherence coefficient e (in log-log scale) for the GDM. The map is the perturbed cat of Eq. l|3j with a = b = 2 and k = 0.01, 
N = 800. The (solid) line is an approximate (e <C 1) analytic calculation F^ = 4 '='^P( j-^^^^^^^_2^2j^gj^-)2^^ ^ji^^ ^ ^ The horizontal (dashed) 

line corresponds to A = ln[(3 + 2\/2)]. The inset shows the same plot in linear scale where the e-independent regime with Fg = A is clearly 
observed. All the results were obtained from the evolution of a single initial state [no average was performed]. 



the behavior of in the case of quantum maps is the same. Af- 
ter a critical value the decay rate saturates to a constant value 
given by A. However, as we show in Fig. |2] [right], for e small 
the dependence of Fg is nowhere near quadratic. For the GDM 
[Fig.|2] right] if £ is very small, of order then the proba- 
bility of applying any translation is negligibly small. Thus for 
e < 1 /N there is no decoherence and the purity remains con- 
stant and equal to unity. For larger decoherence strengths, the 
purity decays exponentially but the dependence of is not 
quadratic. We remark that all the calculations done for the pu- 
rity do not need any kind of averaging. Fig. |2]was obtained 
using a single Gaussian initial state. 

We can derive an approximate analytic expression for the 
small £ regime. If we assume (9,trp^ = Atr(p„+i)^ = tr(p^^j) — 
tr(p2) = -r£tr(p,^), then from Eqs. (0 and I©, if £ < 1, we 
have 

p' « Cs{0,0)p+Ce{0,l)ToApTli+Cs{l,0)Ti,opTlo 

+c, (- 1 , o)r_ 1 ,op ri 1 + (0, - 1 )ro^_ ip roV( 10) 

We want to take the square of the trace, so the first approx- 
imation we take is tr(p7]jp7j'y),-.y=o.i « tr(p^), and we ne- 
glect higher order terms as well as higher order translations 
(even We also take into account the fact 

Ce{q,p) is symmetric around {q = 0, p = 0). Thus we have 
trp'2_trp2 fa (c£(0,0)2 + 4ce(0, 1)^ - l)trp2. Now, neglecting 
also higher order terms in the normalization [remember that 
'Lci.pCe{q,p) = 1], we get 

2 ^ exp(-2;rV(£jV)^) +4exp(-2;rV(£jV)2)^ 2 
^ (l+4exp(-2;rV(eA^)^) + ---)^ 

(11) 



For small £ we can of course neglect the terms coming from pe- 
riodic boundary conditions. This expression reproduces very 
well the results obtained numerically [see Fig. |2] right, solid 
line] in the £ « (9(1 /N) region. 

In order to attain the quadratic dependence of for small 
coupling, observed in continuos Hilbert space, C£{q,p) should 
have tails that decay slower than Gaussian, i.e. long distance 
correlations in phase-space. We can for example take a well 
known decoherence channel for quantum information process- 
ing, the depolarizing channel (DC) |26|, which is also a convex 
sum of unitaries and can be simply written in terms of transla- 
tions in phase space [,23 J 

= (i-£)p + -^ ^ r^^pri-p (12) 

In Fig|3][left] we show the decay rate Fg for the DC. Following 
a similar reasoning as the one followed to obtain Eq. (fTTT l we 
get, for £ < 1, Fg = 2£ [see Fig[3] left, solid line]. The DC is an 
extreme case to consider as phase space decoherence because 
it is highly non-local: with the same probability it implements 
every possible translation T^^ (q,p ^ 0). Therefore, there is no 
reason to expect a Lyapunov regime in this case. In fact for £ 
close to 1, the dynamics is dominated by the environment. The 
absence of Lyapunov regime (or any £ independence) is clearly 
appreciated in Fig [3] [left] . The non-locality of DC has also 
devastating effects on the entangling power of the algorithms 
that implement chaotic maps ||27l . 

To reproduce the FGR quadratic regime we thus need a deco- 
herence model which is peaked at C£(0,0) and which has poly- 
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Fig. 3: Decay rate of the purity as a function decoherence coefficient e (in log-log scale) for decoherence models: (left) DC; (right) LDM. 
The map is the perturbed cat of Eq. ([3) with a = b = 2 and k = 0.01, A' = 800. The solid lines shows the fit for small e. For the DPC (left) 
the solid curve was obtained analyticaly and is F^ =2£. For the LDM [right] the fit is aproximately Fg oc jj^g dashed line in both plots 
corresponds to A = ln[(3 + 2\/2)]. The inset in both panels displays the same plot in linear scale where the absence of an e-independent 
Lyapunov regime is clearly observed for both DC and LDM cases. 



nomially decaying tails. We propose to take a Lorentzian 
1 X eN 

Cs{q,p) = — Y. 7 2 V (13) 

''^-^ ((f ) + ^<i-^jf + {p-m") 

with A the proper normalization for Y,q,pC£{q,p) = 1- We will 
call this case Lorentz decoherence model (LDM). The sum is 
done to account for the periodicity of the torus (theoretically 
.X — > °o, practically x is an integer much larger than 1). Eq. (|7|l 
with C£{q,p) given by Eq. (fTST l defines a random process with 
Lorentzian weight. We can relate this to superdiffusion by 
Levy flights. Long tail decoherence was also considered in 
Ref. ||28l where it was shown that the decoherence rates can 
be tuned to power law decay in cold atom experiments. In 
Fig. [3] [right], we show Fg for the LDM. The quadratic depen- 
dence is clearly observed. As in the DC model the Lyapunov 
regime is not present. Larger e implies longer Lorentzian tales 
which, when periodized sum up to non negligible non-local ef- 
fects all over phase space. This is why for the LDM not only 
is the Lyapunov regime also not present but the decay rate of 
the purity continues to grow indefinitely. To obtain the so- 
called universal behavior - quadratic-FGR growth followed by 
constant-Lyapunov - a very specific model with large tails but 
sufficiently localized is needed. A combination of both GDM 
and LDM, so that the former dominates at larger e and the latter 
dominates for smaller £ would yield both the FGR regime and 
the Lyapunov regime. Decoherence combining both Gaussian 
and Lorentzian processes was studied e.g. in [29 J . 

To summarize, the LE and the purity for systems with fi- 
nite dimensional Hilbert space has been analyzed. We have 
shown that though they can exhibit qualitative similarities, they 
are fundamentally very different: the small coupling regime 
for the purity is not quadratic but depends on the environment 
model. Moreover, while the large perturbation regime for the 



LE can present high amplitude oscillations around the classical 
Lyapunov exponent depending on the LDOS, for the purity it 
depends decidedly on the type of environment. Only environ- 
ments that act locally in phase space exhibit the - independent 
- Lyapunov regime. Thus, we remark that the LE and the purity 
provide intrinsically different information. 
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